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ABSTRACT 

We present [Cn] observations of 20 strongly lensed dusty star forming galaxies at 
2.1 < z < 5.7 using APEX and Herschel. The sources were selected on their 1.4 mm 
flux (Si. 4 mm > 20 mjy) from the South Pole Telescope survey, with far-infrared (FIR) 
luminosities determined from extensive photometric data. The [Cii] line is robustly 
detected in 17 sources, all but one being spectrally resolved. Eleven out of 20 sources 
observed in [Cii] also have low-J CO detections from ATCA. A comparison with mid- 
and high-J CO lines from ALMA reveals consistent [Cii] and CO velocity profiles, 
suggesting that there is little differential lensing between these species. The [Cii], low- 
J CO and FIR data allow us to constrain the properties of the interstellar medium. 
We find [Cii] to CO(1-0) luminosity ratios in the SPT sample of 5200 =b 1800, with 
significantly less scatter than in other samples. This line ratio can be best described 
by a medium of [Cii] and CO emitting gas with a higher [Cii] than CO excitation 
temperature, high CO optical depth tco(i-o) ^ 1? and low to moderate [Cii] optical 
depth T[cii] < 1. The geometric structure of photodissociation regions allows for such 
conditions. 

Key words: galaxies: high-redshift - galaxies: ISM - galaxies: starburst - infrared: 
galaxies - submillimeter: galaxies 
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1 INTRODUCTION 

The discovery of (sub)millimeter-selected dusty star form- 


ing galaxies (DSFGs) at high redshifts (e.g. 

Smail et al. 

1997 Hughes et al. 1998; Barger et al. 

1998 

1 fundamen- 

tally changed our view of galaxy formation and evolution. 


DSFGs are massive (M* ~ lO n M 0 ; e.g. 


Hainline et al. 


2011 

Michalowski et al. 20121 and gas-rich (M ffas ~3-5x 

10“ 

Vt 0 ; e.g. Greve et al.]2005 

Tacconi et al. 2010 

Bothwell 

et al 

2013b 1, and have star for 

nation rates > 1000 M 0 yr 1 
e properties of these galaxies 

(e.g. Chapman et al.|2005). Tli 


remain a challenge for conventional galaxy formation mod¬ 
els (e.g. Baugh et al. ||2005 Lacey et al.||2010 Benson|[2012 


Hayward et al.| 2013). 

Thanks to the availability of new space- and ground- 
based sub-millimetre facilities, our knowledge of the in¬ 
terstellar medium (ISM) in massive gas-rich galaxies at 
high redshift has dramatically improved in the past decade 


(Solomon & Vanden Bout 2005 Carilli & Walter 2013 Casey 


et al.|20l4 1. The most commonly used lines for studying the 
ISM in DSFGs at high redshift are the rotational transitions 
of carbon-monoxide (CO) and the [Cn] A158 /im fine struc¬ 
ture line. The latter is the most important cooling line in 
the ISM ( |Dalgarno~fe M cCray| 1 972] ), and traces neutral gas 
exposed to ultraviolet photons from young stars. The [Cii] 
line can therefore be used to probe the stellar radiation field 
and how it affects the physical conditions of the gas. The 
bulk of the [Cii] emission line (70% in Stacey et al. 1991a 


2010) is believed to originate from photodissociation regions 


(PDRs), and the remainder from X-ray dominated regions 
(XDRs), cosmic ray dominated regions (CRDRs), ionised re¬ 
gions (Hii regions) (Meijerink et al. 20071, low density warm 
gas and/or diffuse HI clouds (|Madden et al.||1997|). 


The [Cii] fine-structure transition (i-'res 1 ? = 

1900.54 GHz) is nearly unobservable from the ground 
at 2 < 1 due to strong atmospheric absorption. The only 


(Crawford et al. 1986 Stacey et al. 

1991a 1 or space based 

observatories (Malhotra et al. 2001 


Brauher et al. 2008 

Diaz-Santos et al. [2013 

|Farrah et al. 

2013; Sargsyan et al. 

2014 De Looze et al. 

2014 1 . Submillimetre atmospheric 


windows provide some access to the line from the ground 
at 2 > 1, with atmospheric transparency continuing to 
improve towards higher redshifts (longer wavelengths). As 
a consequence, the first high -2 [Cii] detection was reported 
for the 2 = 6.42 quasar host galaxy SDSSJ1148+5251 a 


decade ago (Maiolino et al. 2005). The number of [Cii] 
detections has been steadily increasing since then, thanks 
to facilities such as the Caltech Submillimeter Observatory 
(CSO), the Submillimeter Array, the IRAnQ Plateau de 
Bure Interferometer, the Atacama Pathfinder Experiment 
(APEX), the Herschel Space Observatory , the Combined 
Array for Research in Millimeter-wave Astronomy, and the 
Atacama Large Millimeter/submillimeter Array (ALMA). 
At 1 < 2 < 2, 21 [Cii] detections were made using the 
redshift (2) and Early Universe Spectrometer (ZEUS) on 


CSO (Hailey-Dunsheath et al. 2010 Stacey et al. 2010 


Brisbin et al. 20151, while Herchel detected three [Cii] 
lines at 1.5 < 2 < 3 (Ivison et al. 2010 Valtchanov et ah] 


2011 George et al. 2013[ ). Seventeen have been added to 
the number of 2 > 4 [Cii] detections in the past decade 
|Maiolino_et^aL 2005[ |2009 |Iono et al. |2006| Wagg et al. 


2010 Ivison et al. 2010 De Breuck et al. 2011 Cox et al. 


2011 

Swinbank et al. 2012 

Venemans et al. 

2012J Walter 

et al. 

2012 


Riechers et al. 

2013| Wang et al. 

2013; Rawle 

et al. 

2014 


De Breuck et al.||2014j Neri et al. 

2014). Many 


of these objects have been selected based on the presence of 


a luminous active galactic nuclei (AGN; e.g. 

Maiolino et al. 

2005 Stacey et al. 2010 Wang et al. 2013 

1 , while others 


have been selected as starburst galaxies (e.g. 

Ivison et al. 

2010; Hailey-Dunsheath et al. 2010; Stacey et al. 

2010 

Cox et al. 2011 Swinbank et al. 2012; Walter et al. 

2012 

Riechers et al. 

2013 

Brisbin et al. 

2015 

|. This leads to a 


heterogeneous sample of high -2 [Cii] detections, containing 
a mixture of AGN and starburst-dominated systems. The 
heterogeneity of the sample complicates the interpretation 
of trends within it. |Stacey et a l. (2010) suggest, based on 
the [Cii] /FIR and [Cii] to CO (1-0) luminosity ratios, that 
[Cii] emission originates mainly from PDRs and that the 
ISM and stellar radiation field in these 2 ~ 1 — 2 galaxies 
resemble that observed for local starburst systems. These 
studies have followed the conclusion of ICrawford et aid 


(1985 ) that the [Cii] emission is optically thin or reaching 
unity opacity (r < 1); however, this was recently challenged 
by Neri et al. (20141, who argued for optically thick [Cii] . 

Here we present [Cii] observations of 20 gravitationally 
lensed DSFGs in the redshift range 2 ~ 2.1 — 5.7 discovered 


by the South Pole Telescope (SPT; Carlstrom et al. 2011 


Vieira et al. 2010). These 20 sources are a subset of those 
selected from the first 1300 deg 2 of the 2500 deg 2 SPT-SZ 
survey. Followup observations with ALMA have provided 
spectroscopic redshifts for all of these objects (|Weifi et al. 


[20l3l [YleTra et al.||2013| ). We also include two new sources 
observed in Cycle 1 (see Appendix[A]for details). The bright 
(*Si. 4 mm > 20mJy) flux selection of the SPT sample ensures 
that virtually all sources will be gravitationally magnified, 


with a bias towards 2 > 2 (Hezaveh & Holder 2011 Weifi 
et al.|[2013). The magnified emission allows us to study the 


ISM in these DSFGs in greater detail, using fine structure 
and molecular lines such as [Cii] and CO. By including low- 
J CO observations for 11 sources in our analysis, we deter¬ 
mine the physical state of the ISM by studying the [Cii] and 
CO(l-O) line intensity ratios. 

This paper is organised as follows: in [J2]we describe the 
[Cii] and CO observations, and the results for these observa¬ 
tions are given in m In EJwe present our analysis, and dis¬ 
cuss the implications in 95] Our conclusions and summary 
are given in E Throughout this paper we adopt the cos¬ 
mology: Ho = 71km s -1 Mpc -1 , Da = 0.73 and D m = 0.27 
(Komatsu et al.|2011i. 


Institute for Radio Astronomy in the Millimeter Range 


2 OBSERVATIONS 

2.1 Supporting ALMA and ATCA observations 

The 20 DSFGs presented here are a subset of 100 strongly 
lensed DSFGs selected over the 2500 deg 2 SPT-SZ survey. 
See Table 1 in Weifi et al. ( |2013| ) and Table Al in Ap- 
pendix[X]for the full names and positions. Mid- and high-J 
CO rotational lines (i.e. CO(3-2), CO(4-3), CO(5-4) and/or 
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CO(6-5)) were detected for 23[2]DSFGs with the Atacama 
Large Millimeter/submillimeter Array (ALMA). The red- 
shifts of all our [Cn] targets are determined by one or more 
CO lines plus the [Cn] line itself and are therefore robust. 

In addition, low-J CO emission lines (CO(l-O) or 
CO(2-l)) were observed with the Australia Telescope Com¬ 
pact Array (ATCA) for 11 of the SPT DSFGs for which 
[Cn] observations are presented in this paper 
[2013] Aravena et al. in prep). Absolute flux 
the ATCA data is estimated to be accurate to within 15%. 
Details of the observations, fluxes and associated uncertain¬ 
ties will be presented in a forthcoming paper by Aravena et 
al. 


Aravena et al. 
calibration of 


2.2 APEX/FLASH 


We targeted all galaxies in the SPT DSFG sample with 
known redshifts that place the [Cn] line at frequencies 
which are observable with good atmospheric transparency 
using the First Light APEX Submillimetre Heterodyne re¬ 


ceiver (FLASH, Heyminck et al. 2006). Eleven sources at 

4.2 < 2 < 5.7 were observed in the 345 GHz channel 
between 2012 August and 2014 June, and six sources at 
3.1 < 2 < 3.8 were observed with the 460 GHz channel be¬ 
tween 2013 March and August during Max Planck time. All 
observations were done in good weather conditions with an 
average precipitable water vapour < 1.0 mm, yielding typi¬ 
cal system temperatures of 230 and 170 K for the 345 and 
460 GHz observations, respectively. The beam sizes/antenna 
gains are 22.0”/40 Jy K^ 1 and 13.5”/48 Jy K _1 for the low¬ 
est and highest observed frequencies of the [Cn] line, re¬ 
spectively. The beam size is much larger than the observed 
Einstein radii of these sources and thus they are unresolved 


(Vieira et al. 2013 Hezaveh et al. 2013). The 82 hours of 


observations were done in wobbler switching mode, with 
switching frequency of 1.5 Hz and a wobbler throw of 50” 
in azimuth. Pointing was checked frequently and was found 
to be stable to within 2.5”. Calibration was done every 
~ 10 min using the standard hot/cold-load absorber mea¬ 
surements. The data were recorded with the MPIfR Fast 
Fourier Transform spectrometers (FFTS; Klein et alj2006) 
providing 4 x 2.5 GHz of bandwidth to cover the full 4 GHz 
bandwidth in each of the upper and lower sidebands of the 
sideband-separating FLASH receiver. 

The data were processed with the Continuum and Line 
Analysis Single-dish Software (CLASfQ). We visually in¬ 
spected the individual scans and omitted scans with unsta¬ 
ble baselines, resulting in < 10% data loss. We subtracted 
linear baselines from the individual spectra in each of the 
two FFTS units, and regridded to a velocity resolution of 
~ 90kms _1 in the averaged spectra. On-source integration 
times were between 1.5 and 5 hours. Table [I] summarises the 
line intensities, and Figure [l] shows the spectra. 

We detect [Cii] emission in 16 out of 17 sources ob¬ 
served with FLASH. The only non-detection is in the high¬ 
est redshift source, SPT0243-49 at 2 = 5.699. This source 
has an unambiguous redshift confirmed by two high-J CO 


lines with ALMA ( Weifi et al.|20 f31 and a CO(2-l) line with 
ATCA (Aravena et al. in prep.). The observed L[cii]/Lfir 
ratio (< 1.3 x 10 -3 ) is close to the median of our sample (see 
Table [I]), suggesting that a [Cii] detection for this source is 
feasible with a moderately deeper integration. 


2.3 Herschel /SPIRE 

For three 2 ~ 2 sources in the SPT sample (see Table [I]), 
the [Cii] line falls in a frequency range (500-610 GHz) 
where the atmosphere is opaque. We thus observed these 
source s with the SPIRE Fourier Transform Spectrometer 
(FTS; Griffin et al.|2010 1 onboard Hersche Q ( |Pilbratt et al. 
20101. For more detail about the observations and results for 
SPT0538-50, see Bothwell et al.| ( |2013a| ). The observations 
of SPT0551-50 and SPT0512-59 were carried out on 2013 
March 2 in single-pointing mode using both the short wave¬ 
length (SSW) and long wavelength (SLW) bands covering 
194-313 (im and 303-671/rm. The observations were done 
in high spectral resolution mode equivalent to 0.04 cm -1 
(1.2 GHz) with 100 repetitions, resulting in an on-source in¬ 
tegration time of 13752 s (3.8 h) per source. 

The data were reduced and calibrated using an updated 
SPIRE FTS pipeline in the Herschel reduction tool HIPE 
vll, which includes all detectors in the observation and uses 
new calibration files. The detection of lines fainter than 1 Jy 
is very challenging because thermal emission from the warm 
optics contributes as much as 1000 Jy at 1000 GHz. We sub¬ 
tract the average of the off-target pixels in order to remove 
this excess emission from the telescope. Another possibility 
is to subtract a ‘dark sky’ observation made on the same 
day with the same exposure time, but in our case the noise 
level was lowest by subtracting the average of the off-target 
pixels. As the model of the telescope has an uncertainty 
of 0.1%, there remains a residual continuum uncertainty of 
~ 1 Jy in the continuum. We therefore subtract an addi¬ 
tional second order polynomial from the SSW and the SLW 
spectral part separately, and look for the [Cii] line at the 
expected frequency (Valtchan ov et al.|2011| ). Figure [2] shows 
the resulting spectra. 


3 RESULTS 

Figures [I] and 2 show the [Cn] emission lines and Table [l] 
lists the [Cii] luminosities obtained from APEX/FLASH and 
Herschel SPIRE FTS observations, along with far-IR lumi¬ 
nosities (Lfir). We have 17 [Cn] emission line detections 
(16 with FLASH and one with Herschel SPIRE FTS) and 
three non-detections (one with FLASH and two with Her¬ 
schel SPIRE FTS). 


3.1 Velocity profiles and line fluxes 

Despite an increasing number of high -2 [Cn] detections in 
the literature, only a few of these have sufficient spectral res¬ 
olution and S/N and supporting data to compare line pro- 


One source SPT0538-50 ( jGreve et al.|[2012| |BothweIl et al.] 
, was not observed with ALMA, 
http: //uw. iram. f r/IRAMFR/GILDAS/ 


4 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 
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files with other bright lines like CO (e.g. Rawle et al. |2014 1. 
This is the first sample of sources with spectrally resolved 
data with S/N^ 3 in both [Cn] and CO making it possible to 
compare the shapes of velocity profiles. In the following, we 
consider only the 17 sources observed with FLASH, as the 
3 sources observed with SPIRE are not spectrally resolved. 

Figure [l] shows the velocity profiles of the [Cii] lines 
compared with the mid- and high- J CO lines observed with 
ALMA. The CO lines have been scaled to match the peak 
flux of the [Cii] line, in order to facilitate the comparison 
of the velocity profiles. We first fit the CO and the [Cii] 
lines independently with single Gaussian functions. We ac¬ 
cept the single Gaussian fit if the reduced y * 2 does not exceed 
1+5 x V^ridof, where ridof is the number of degrees of free¬ 
dom. This quantity corresponds to 5 standard deviations of 
the \ 2 distribution. If the single Gaussian fit does not match 
the above criterion, we use a double Gaussian function, i.e. 
two Gaussian functions with displaced central positions. The 
double Gaussian function is sufficient to describe the line 
profiles that do not match a single Gaussian. In practice, 
this happens in SPT0103-45 and SPT0418-47, which have 
lines that display a slight asymmetry on the red side of the 
lines; the velocity difference between the CO and [Cii] peaks 
is <150kms -1 . 

We then simultaneously fit the single or double Gaus¬ 
sian profile (depending on what is necessary to fit the profiles 
individually) to the CO and [Cii] velocity profiles assuming 
the profiles are similar with just one free scaling parameter. 
This allows us to test if the two profiles are consistent or not. 
Only the [Cii] and CO lines for SPT0532-50 have different 
velocity profiles. The remaining 16 sources have consistent 
line profiles with \ 2 4 * * * < 1 + 5 x *j2n& 0 {. 

We obtain the line widths (FWHM) listed in Table [l] 
by fitting a single Gaussian. To test the reliability of our 
method we also fit the spectra by taking velocity-weighted 
moments, and find fully consistent results. The resulting 
FWHMs are in the range ~ 210 — 820km/s. This is compa¬ 
rable to the typical CO line width of ~ 460kins -1 found for 
SPT DSFGs by stacking 22 spectra (Spilker et al. 2014). 
Nine out of 20 sources have FWHM > 500km/s, which 
is large compared to other high-z [Cii] detections (e.g. 
360km/s on average in the sample of |Wang et al.|2 013). 

We obtain the velocity-integrated fluxes for the 
FLASH/APEX detections by summing the observed line 
profiles over the 3<r limits obtained from the single Gaus¬ 
sian fits. The [Cii] apparent luminosities range from 1.4 to 
9.2 x 10 10 Lq which is 1-2 orders of magnitude higher than 
their mid- and high-J CO luminosities (~ 3 — 30 x 10 8 * Lq), 
where both are uncorrected for lensing. 

The [Cii] detection in SPT0512-59 with SPIRE FTS 
(see Figure [2] right) confirms the redshift at z = 2.234. We 
determine the integrated line flux in the same manner as 


Valtchanov et al. (20111 by fitting the emission lines with a 


sinc-function, and calculate the RMS using the standard de¬ 
viation of each channel within ±5000 kms -1 of the expected 
line centroid (see Valtchanov et al. 2011 for more details). 
Our SPIRE FTS observation for SPTO551-5C0does not 


5 Note that is is one of the rare cluster lenses in the SPT sample 
| [Vieira et akpObjl). This, however, should not have any effect on 
the detectability of [Cii] in this source. 


detect any [Cii] emission (see Figure [ 5 ] left). The redshift 
was reconfirmed by weak CO (1-0) emission observed with 
ATCA (Aravena et al. in prep.). We therefore take the 30- 
upper limit of the expected line peak to be 3 times the RMS 
noise. 


3.2 Lensing 


The lensing magnification of the SPT DSFGs allows us to 
study the ISM in galaxies at high redshifts, but also intro¬ 
duces the possibility of differential lensing. The compact¬ 
ness and location of a region relative to the lensing caustic 
determines the magnification of the emission. Differential 
lensing amplification may thus occur between compact and 
extended emitting regions, or components occupying differ¬ 
ent regions (e.g. Hezaveh et al.|2012 l. From observations of 
the Milky Way and nearby galaxies, we might expect the 
[Cii] emission to originate from more extended and diffuse 
media than the more optically thick low-J CO emission in 
giant molecular clouds (GMCs) (Fixsen et ah]|1999 1. If the 
[Cii] was dominated by emission from such diffuse regions, 
it could be subject to differential lensing compared to the 
more compact CO emission. However, the similar [Cii] and 
CO velocity profiles suggest that such differential lensing is 
not significant. 

Modelling the lensing magnification factor is of great 
importance, and has been performed for some of the DSFGs 
in this sample using ALMA data. The lens modelling is per¬ 
formed in the (u, v )—plane to properly represent the param¬ 
eter uncertainties in the interferometric data (see |Hezaveh 
et al.| 2013 for more details). The four sources with lens mod¬ 
els have lensing magnification factors p between ~ 5.4 — 21.0 
( |Hezaveh et al.|[2013[ ) and a mean of (p) = 14.1. In cases 

where the lensing magnification factor is unknown, we use 

the mean magnification factor to estimate the intrinsic lu¬ 

minosity. We then conservatively choose to span the uncer¬ 
tainty on the lensing magnification factor from the smallest 
(p = 5.4) to the largest (p = 21.0) giving the mean mag¬ 
nification an uncertainty of 7.8. We assume that the mag¬ 
nification factors derived from the dust continuum are also 
appropriate for the line emitting gas. A lens modelling anal¬ 
ysis of the remainder of the ALMA Cycle 0 imaging data is 
currently under way. 


4 ANALYSIS 

4.1 Comparison sample of nearby and distant 
galaxies 


We compare the SPT sample with that of Gracia-Carpio 

et al. (in prep), which is comprised of 333 sources. These 

sources include LIRGs and ULIRGs from the Great Obser¬ 


vatories All-sky Survey (GOALS, Dfaz-Santos et al. 20131 


and normal and Seyfert galaxies from |Brauher et al.| ( 20081 

with spatially unresolved [Cii] detections integrated over the 
entire individual galaxies. The sample contains 308 sources 
at z < 0.4 and 25 at z > 1. In addition, we have searched 
the literature for additional z > 1 DSFGs with [Cii] obser¬ 
vations, which we list in Appendix [B] In constructing this 
low and high-s comparison sample, we have paid particular 
attention to ensure that the photometric data are integrated 
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Figure 1. Comparison of the velocity profiles of the [Cn] lines detected with APEX/FLASH and SPIRE FTS (grey filled profiles) and 
mid-J CO line observed with ALMA in Cycle 0 (coloured lines; [Weifi et al.|2013[ ).The CO lines have been scaled to match the [Cnjpeak 
flux. The similarities between the CO and [Cn]lines for individual sources suggest that the spatial distributions are similar and differential 
lensing is not significant. 
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Figure 2. Herschel SPIRE FTS spectra for SPT0551-50 and SPT0512-59. Left: Non-detection of the [Cn] emission line for SPT0551-50. 
Right: Detection of [Cn] emission line for SPT0512-59. The continuous red curve is the sinc-function used to fit the [Cn] line (see 
|Valtchanov et al.|2011 ). 


Source) 

z 

SdV { cm 

[Jykm/s] 

dV 

(FWHM) 

[km/s] 

L{cii]/10 10 

[K km/s/pc 2 ] 

T[cii] 

[10 10 L 0 ] 

Tfir 

[10 13 L 0 ] 

T d 

[K] 

L [CII] 

^FIR 

[io- 4 ] 

Instrument 

Time 

[h] 

SPT0551-50 1 

2.123 

< 180 (3a-) 

— 

< 13.4 

< 3.0 

1.1 ± 0.1 

27.2 ± 1.0 

< 26.6 

SPIRE FTS 

3.8 

SPT0512-59 

2.234 

227 ± 43 

— 

18.4 ±3.5 

4.0 ±0.8 

2.8 ±0.2 

33.2 ± 1.2 

14.4 ±2.9 

SPIRE FTS 

3.8 

SPT0538-50 1 

2.782 

< 465 (3<r) 

— 

< 81.9 

< 18.0 

5.8 ±0.3 

36.9 ± 1.4 

< 31.0 

SPIRE FTS* 

3.8 

SPT0103-45 

3.090 

125 ± 17 

304 ± 47 

17.5 ±2.3 

3.8 ±0.5 

3.4 ±0.2 

33.5 ± 1.1 

11.3 ±1.6 

FLASH 

2.7 

SPT0550-53 

3.129 

129 ± 25 

719 ± 124 

18.5 ±3.6 

4.1 ± 0.8 

1.6 ±0.1 

34.5 ± 1.8 

25.4 ±5.2 

FLASH 

17.3 

SPT0529-54 

3.369 

217 ± 18 

823 ± 92 

35.1 ± 3.0 

7.7 ± 0.7 

3.0 ±0.2 

33.1 ± 1.2 

25.7 ± 2.8 

FLASH 

3.8 

SPT0532-50 

3.399 

113 ± 18 

767 ± 124 

18.6 ±3.0 

4.1 ± 0.7 

6.5 ± 0.4 

37.9 ± 1.4 

6.3 ± 1.1 

FLASH 

4.6 

SPT0300-46 

3.596 

41.5 ±10.4 

583 ± 138 

7.5 ± 1.9 

1.6 ±0.4 

3.3 ±0.2 

39.2 ± 1.5 

5.0 ± 1.3 

FLASH 

12.8 

SPT2147-50 2 

3.761 

80.5 ± 11.7 

329 ± 56 

15.5 ±2.3 

3.4 ±0.5 

3.2 ± 0.2 

41.4 ± 1.7 

10.7 ± 1.7 

FLASH 

5.0 

SPT0418-47 2 

4.224 

127 ± 10 

347 ± 29 

29.5 ± 2.3 

6.5 ± 0.5 

5.9 ±0.4 

47.3 ± 2.5 

11.0 ± 1.1 

FLASH 

1.5 

SPT0113-46 2 

4.232 

91 ± 19 

619 ± 132 

21.1 ± 4.4 

4.6 ± 1.0 

2.1 ± 0.1 

32.9 ± 1.4 

22.1 ± 4.8 

FLASH 

4.5 

SPT2311-54 

4.281 

46.4 ± 5.3 

360 ± 44 

11.0 ± 1.2 

2.4 ±0.3 

3.3 ±0.3 

43.3 ±3.3 

7.3 ± 1.1 

FLASH 

3.0 

SPT0345-47 2 

4.296 

63.7 ± 8.3 

810 ± 200 

15.2 ±2.0 

3.3 ±0.4 

9.2 ± 0.8 

51.8 ±3.2 

3.6 ±0.6 

FLASH 

2.3 

SPT2103-60 2 

4.435 

129 ± 18 

780 ± 125 

32.2 ±4.4 

7.1 ± 1.0 

3.4 ±0.2 

39.2 ± 1.5 

20.8 ±3.1 

FLASH 

1.5 

SPT0441-46 2 

4.477 

42.5 ± 10.6 

581 ± 162 

10.8 ±2.7 

2.4 ±0.6 

3.7 ± 0.2 

39.9 ± 1.9 

6.6 ± 1.7 

FLASH 

3.2 

SPT2146-55 2 

4.567 

39.0 ± 9.0 

302 ± 62 

10.2 ±2.4 

2.2 ± 0.5 

2.7 ±0.3 

39.2 ±2.0 

8.3 ± 2.1 

FLASH 

3.1 

SPT2132-58 

4.768 

34.9 ± 6.9 

212 ±43 

9.7 ± 1.9 

2.1 ± 0.4 

3.1 ± 0.3 

39.5 ± 1.9 

6.9 ± 1.5 

FLASH 

2.1 

SPT2319-55 

5.293 

19.1 ± 3.2 

198 ± 34 

6.3 ± 1.1 

1.4 ±0.2 

2.5 ± 0.2 

42.0 ±3.1 

5.4 ± 1.1 

FLASH 

9.5 

SPT0346-52 2 

5.656 

63.3 ± 8.7 

502 ± 72 

22.8 ±3.2 

5.0 ±0.7 

12.3 ±0.5 

52.4 ± 2.2 

4.1 ± 0.6 

FLASH 

1.4 

SPT0243-49 2 

5.699 

< 51 (3a-) 

— 

< 21.0 

< 4.5 

3.3 ±0.3 

35.3 ± 1.6 

< 13.6 

FLASH 

2.2 


Table 1. Observed [Cii] and FIR properties. All luminosities are uncorrected for the lensing amplification. The upper limits for the 
velocity integrated fluxes given for the non-detections are obtained by assuming the FWHM observed for the CO lines. The total 
integration time with SPIRE FTS and FLASH is 92 hours. The integration time per sources is given in the last column, 
f Full source names are listed in Table 1 of |Weifl et al.||2Q13| or in Table [Al] 

1 Has CO(1-0) observations (Aravena et al., in prep.). 

2 Has CO(2-1) observations (Aravena et al., in prep.). 

★ See |Bothwell et al.| ( |2013a| ) for more details. 


over the full galaxies, as all line and continuum data for the 
SPT sample are also integrated values. 


4.1.1 Conversion to CO(1-0) luminosities 

Like the SPT sample (§ ED- many of the sources in the 
comparison sample also have published CO observations. In 
cases where the CO (1-0) emission lines have not been ob¬ 


served, Gracia-Carpio et al. (in prep) derive the Lco(i-o) by 
converting the observed mid- J CO luminosities to L co(i-o) 
using scaling factors from Stacey et al. (2010). When mul¬ 
tiple J > 1 CO lines were observed, we take the average 
of the scaled Tco(i-o)- These scaling factors are based on 
previous studies involving several rotational lines from both 
nearby galaxies, ULIRGs and high -2 galaxies (|Stacey eF al. 


2010), which allow for the assumption of fixed CO line ra- 
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tips up to CO(4-3)/ CO(l-Q). The conversion factors used by 
Stacey et al. (2010) assume an integrated line flux (Wm -2 ) 


ratio of CO(2- 1)/CO(1-0) = 7.2, equivalent to 90% of the 
thermalised optically thick emission (i.e. a brightness tem¬ 
perature ratio of 0.9), CO(3-2)/CO(2-l) = 3.0 (90% of the 
thermalised optically thick emission) and CO(4-3)/CO(2-l) 
= 6.4 (80% of the thermalised optically thick emission). 

The assumed CO(2-l) to CO(l-O) brightness temper¬ 
ature ratio of 0.9 for DSFGs is in agreement with observa¬ 
tions of [Bothw ell et al.| f2C)13b), who find a ratio of 0.85 for 
32 luminous submm galaxies at z ~ 1.2 — 4.1. The [CO(2- 
1)/CO(1-0)] ratio has also been observed for four normal 
star forming galaxies by |Aravena et a l.j ( |2014[ ), resulting in 
a slightly lower [CO(2-1)/CO(1-0)] brightness temperature 
ratio of 0.7±0.16 for four BzK galaxies at 2 ~ 1.5 —2.2. Even 
though the ratio is consistent with Bothwell et al.| (2013b) 
within the uncertainties, a lower [CO(2-1)/CO(1-0)] ratio is 
expected for normal star forming galaxies than for starburst 
galaxies, because the molecular gas in starburst galaxies is 
expected to be more highly excited. Spilker et al. (2014) find 
an average CO(2-1)/CO(1-0) brightness temperature ratio 
of 1.1 ± 0.1 for 22 of the SPT DSFGs by stacking ATCA 
spectra, after scaling them by their 1.4 mm continuum flux 
density. Given the range of values of these three methods, 
and for consistency with previous literature, we adopt the 
scaling factor of 0.9 from Stacey et al. (2010). 


4-1.2 FIR luminosities 


We obtain the FIR luminosities Lfir and the dust temper¬ 
ature Th for each source in the SPT DSFG sample by fit¬ 
ting the well-sampled spectral energy distributions (SEDs) 
with a greybody law fixing the emissivity index (/3) at 2.0, 
Uo = 100 and fitting A res t < 50/rm, following |Greve et al.| 
12012) and Strandet et al. (in prep). We integrate the SED 
between Arest = 42 — 500/rm in the rest frame to obtain 
Lfir. Our SPT DSFGs all have 7 photometric data points 
covering observed wavelengths from 250 to 3000 /im (e.g. 
Weifi et al.||2013l. This allows for a uniform determination 
of the FIR luminosity using a parametrised SED fitting] 
The SPT sources, which have well-sampled SEDs, thus have 
a smaller uncertainty in Lfir than sources with poor pho¬ 
tometric coverage (see Figure [3]. 

To compare the low -2 and high -2 sample with the SPT 
sources in a consistent way, we compiled published FIR pho¬ 
tometry (Gracia-Carpio et al. in prep; Appendix E , and 
derived Lfir and Tj using the same fitting code we used 


for the SPT sources (Greve et al. 2012 1 . IRAS and ISO 


data are available for sources published by Brauhcr et al.| 
([2008]), and FIR data are also available for a large number of 
the GOALS sources at NASA/IPAC Extragalactic Database 
(NEE0. The comparison sample contains 165 sources with 
sufficiently good photometry to derive Lfir and T& using 


the method in Greve et al. (2012); 14 of these are 2 > 1 
sources (see Figure |3j). Sources with insufficient available 


6 Other studies (e.g. jStacey et al.|2010||Helou et ;i 1 . 11 988[ i use the 
two band definition for Lfir: Lfir = 1.26 x 10“ 14 x (2*58/60 + 
/ioo)[W m — “] which is equivalent to the 42.5-122.5 pm luminosity 
l Hclou et al. 1988 l. 

' http:/7ned.ipac.caltech.edu 



Z 


Figure 3. Observed Lfir vs redshift for the 20 SPT sources and 
the comparison sample. No Lfir have been corrected for lensing 
magnification factors. The distribution shows the Malmquist bias 
where high -2 galaxies require either lensing magnification or very 
high intrinsic FIR luminosity of Lfir (7 10 12 Lq to be included 
in the parent sample. The evolution of the luminosity function 
and the smaller comoving volume at low redshifts imply that ob¬ 
jects with similar high intrinsic Lfir are missing from the low -2 
sample. However, the most highly lensed DSFGs may have similar 
intrinsic LpiR,than the most luminous sources in the local sample. 


photometry to derive Lfir using our procedure are not in¬ 
cluded in analyses requiring Lfir but are still included in 
the L[cn]/Lco(i-o) ratios. 

Figure [3] presents the observed Lfir distribution of our 
combined sample. There is a clear gap between the low and 
high -2 samples due to the limited sensitivity of the nearby 
samples, and the atmospheric transparency prohibiting [Cii] 
observations at 2 < 1 from the ground. The selection func¬ 
tion of the high -2 sample is quite complex because many 
sources have been pre-selected to have a good chance of de¬ 
tection in [Cii] (e.g. by having strong dust continuum and/or 
CO emission). In addition, the improved atmospheric trans¬ 
parency at lower frequencies mostly compensates for the dis¬ 
tance dimming in the more distant sources. The evolution 
of the luminosity function and the smaller co-moving vol¬ 
ume at low redshifts imply that objects with similar high 
intrinsic Lfir are missing from the low -2 sample. 

The SPT sources in this sample have FIR luminosi¬ 
ties in the range Lfir= (1.2 — 11.8) x 10 12 Lq. After de- 
magnification the range is Lfir//t = (1.1 — 21.9) x 10 12 Lq, 
with a mean de-magnified Lfir of 5.2 x 10 12 Lq, which is 
similar to the most luminous sources in the low -2 sample. 

The uniform sensitivity of the ALMA mm spectroscopy, 
continuum, low-J CO and [Cii] observations (all with 
S/N>3) lead to a high completeness of the SPT sample. 
This, combined with the lensing amplification factor which 
is about one order of magnitude ( Hezaveh et al.|2013 1, make 
the SPT DSFG sample one of the most representative sam¬ 
ples to date for massive, IR-luminous starburst galaxies at 
high redshifts. 
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4-1.3 AGN content 

The comparison sample contains both starburst galaxies 
and luminous AGN. In the high-redshift sample, we dis¬ 
tinguish the AGN-dominated sources from those without 
any known AGN. These high-redshift AGN are quite rare 
sources, mostly selected based on their bright optical emis¬ 
sion lines (e.g. |Wang et al7||2013[ ) . The SPT DSFGs are se¬ 
lected solely on their lensed 1.4 mm continuum flux, and di¬ 
rect mm imaging and spectroscopy avoids any radio/optical 
identification steps that may introduce biases towards AGN- 
dominated systems. Optical spectroscopy of the SPT DSFG 
sample to derive the redshifts of the foreground lensing 
galaxies has not shown any indications of type-1 or type- 
2 AGN. Strongly obscured type-2 AGN may still be present 
in some SPT DSFGs. Supported by the discussion in j j4.4[ 
we will assume in the following that the AGN contributions 
in the SPT DSFGs are negligible. 


4.2 Observed [Cn] to FIR ratios 

Figure [4] presents the T[cii]/Afir ratio against Lfir for the 
SPT sources and the comparison sample. The typical error 
bar for the low -2 sources (in this and the following plots) 
is illustrated by the cross in the bottom left. This typical 
error bar includes the quoted uncertainties of the lines (e.g. 


Young et al. 1995 Negishi et al. 20011, the absolute and 


statistical uncertainty of the FIR photometry, and of the 
FIR luminosity determined by our own SED fitting. The 
SPT DSFGs have been corrected for lensing amplification 
either using the known lens model (Hezaveh et al. 20131, or 
assuming a mean magnification factor (p) = 14.1 ± 7.8 for 


the sources without a lens model (see section 3.21. 

Sources with Lfir < 10 11 L 0 appear to have a roughly 
constant L[ CII ]/Lfir ratio (~4x 10~ 3 ). At Lfir > 10 11 L 0 , 
the ratio drops to 6 x 10~ 4 for the 2 < 1 sources in the com¬ 
parison sample, similar to what has been reported in pre¬ 
vious studies (e.g. Maiolino et al.|[2009| Stacey et al.|2010 |. 
Whether this is an intrinsic or observational effect due to 
limited sensitivity of the [Cn] and/or FIR photometry in 
the low -2 samples is not clear, and investigating this is be¬ 
yond the scope of this paper. The 2 > 1 sources from the 
comparison sample are scattered over two orders of magni¬ 
tude, which may be due to the heterogeneous mix of the par¬ 
ent samples. The highly complete (82% detections) and uni¬ 
formly observed SPT sources have a smaller scatter and an 
average Ljciij/Lfir ratio of ~ 10 -3 . A Kolmogorov-Smirnov 
test results in a probability of 0.7 that the Ljciij/Lfir val¬ 
ues for Lfir > 10 11 L 0 from the SPT and the low -2 samples 
are drawn from the same distribution. 

Normalising the FIR luminosity by the molecular gas 
mass reduces th e scatter seen in Figure [4] for the Lfir > 


10 11 L 0 sources (Gracia-Carpio et al.|2011). The Lfir/ALh 2 


ratio is expected to be proportional to the star formation effi¬ 
ciency, specifically to the number of stars formed in a galaxy 
per unit molecular gas mass. This ratio has the additional 
advantage that the lensing magnification factor cancels out. 
We estimate the molecular gas masses for the 11 sources for 
which we also have low-J CO line detections (Aravena et al, 
in prep.). We determine the molecular gas mass by: 


A-/h 2 = acoLco(i-o)! 


(1) 



10 


10 


10 
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Figure 4. L[cii]/Lfir vs Lfir f° r SPT sources and the com¬ 
parison sample. As reported in previous versions of this plot, 
the LrQiji/LpiR is anti-correlated with Lrir. In particular at 
L > 10 11 L 0 , the l.\c\ 1] / L|, | ijral.io drops and has a larger spread. 
For the SPT sources without known lensing models, we assume 
a lensing magnification factor of 14.1 (and an uncertainty which 
encompasses the range of 5 to 22 from the known models). The 
typical error bar for the literature sources is represented by the 
black cross in the lower left. The histogram on the right shows 
the distribution of galaxies with L F ir > 10 11 L 0 . 


where aco is the CO-to-H 2 conversion factor. To be con- 


sistent with Gracia-Carpio et al. (2011 

1 , we assume aco = 

0.8M 0 (Kkm/spc^) -1 determined by 

Downes & Solomon 


(1998) for the SPT DSFGs and the low- and high -2 com¬ 
parison samples. Figure [5] plots L[cii]/Lfir as a function of 
Lfir/ALh 2 and shows that the DSFGs lie among the local 
LIRGs. Low -2 sources with Lfir > 10 11 L 0 and the liigh -2 
sources have similar L[chi/Lfir ratios. 


4.3 Dust temperatures 


Malhotra et al. 


(19971, the L[ CII ]/L F ir ra- 


As first shown by 
tio shows a strong anti-correlation with Td. They attributed 
this trend to an increase in the Go /n ratio (far-UV ionising 
field over density) in the hotter, more active galaxies, and 
hence a lower efficiency of gas heating reducing the [Cn] flux 
while increasing the dust temperature. We here revisit this 
trend using the uniformly-derived set of Lfir and Td values 
for the SPT sources and comparison samples, derived from 
our SED modelling (§ 4.11. The top panel of Figure [h] com¬ 
pares the L[cii]/L F ir ratio to Td, with both quantities being 
independent of the lensing magnification. The strong anti¬ 
correlation of these quantities is visible in all three samples. 
However, if the dust and [Cn] emission are coming from the 
same regions, and L[cni is less dependent on Td, a signifi¬ 
cant part of this anti-correlation can be explained with the 
Stefan-Boltzmann law which predicts Lfir oc Td 4 . 

In order to look for a residual correlation, we cancel out 
the Lfir oc Td 4 dependence by plotting L[cu] x Ti 4 /Lfir 
against Td (see bottom panel of Figure [6|. We note that a 
small systematic offset in Td would get propagated as Td 4 . 
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Figure 5. T[cii]/Tfir vs the FIR luminosity normalised by the 
molecular gas mass. The molecular gas mass is derived assum¬ 
ing a conversion factor of «co = 0.8 Mq (Kkms -1 pc" 2 ) —4 . The 
T F ir/A1h 2 ratio is expected to be proportional to the num¬ 
ber of stars formed in a galaxy per unit molecular gas mass 
||Gracia-Carpio et al.|2011| . This Mg 2 normalisation removes the 
uncertainty due to the unknown lensing magnification factors, 
and reduces the scatter seen in Figure [4] but the deficit in the 
L[cii]/Tfir ratio still persists. The typical error bar is repre¬ 
sented by the black cross. 



To test the presence of a correlation taking into account 
the [Cii] upper limits, we use the generalised Kendall’s tau 
method (Lavalley et al. 19921. We find a probability that 
both variables are not correlated of 0.036, 0.602 and 0.151 
for the low- 2 , high -2 and SPT samples, respectively. The 
small displacement of the SPT sources relative to the com¬ 
parison sample in Figure [6] is likely due to Malmquist bias 
and evolution effects resulting in the most luminous sources 
being absent in the low -2 comparison sample (see § 4.1.21. 
While there is marginal evidence for a small positive resid¬ 
ual correlation (especially in the low -2 sample), most of the 
correlation seen in the top panel has canceled out. This con¬ 
firms our assumptions that the dust and [Cii] emission are 
originating from the same regions and that 75 [cii] is less de¬ 
pendent on Td- We therefore conclude that the observed 
correlation is mostly dominated by the Stefan-Boltzmann 
law. 


4.4 Observed [Cii] to CO ratios 

One of the strengths of the SPT DSFG sample is that more 
than half of the sources have both [Cii] and low-J CO de¬ 
tections. Figure [7] plots the [Cii] vs CO(l-O) luminosities for 
these 11 sources, along with high -2 sources in the comparison 
sample^] All high -2 sources fall close to the T[cii]/Tco(i-o) 
~ 4400 relation for local galaxies reported byjCrawford et ah] 


8 Note that Figure ^ presents a larger number of high -2 sources 
than Figures [5] [6] and [8] as it does not involve quantities de¬ 
rived from continuum photometry. The calculation of Lfir and 
Td require photometric data that are unavailable for several of 
the high -2 comparison objects (see section 4.1 and Appendix [ b|. 





T dust [K] 


Figure 6. Top panel: Tjqjjj/Tfir vs Td for the SPT DSFGs 
and the low and high -2 comparison sample. The anti-correlation 
between the T[cii]/Tfir ratios and the dust temperatures is 
seen for both low- and high -2 sources, and is expected because 
the Stefan-Boltzmann law predicts TpiR Td 4 . Bottom panel: 
[Cii] xTd 4 /I/FiRVs Td for the SPT DSFGs and the low and high- 
2 comparison sample. Multiplying the T[cii]/T F ir ratio with Td 4 
cancels out the temperature dependence of the Stefan-Boltzmann 
law. All lensing magnification factors and beam filling factors can¬ 
cel in both panels, and the typical error bar is represented by the 
black cross. 


(|1985 ); Wolfire et al. (19891; Stacey et al. (1991a 20101; 


Swinbank et al. (20121; Neri et al. (20141. This relation 


has been explained in terms of PDR models (Wolfire et al. 
1989| [l993| |Stacey et al.||1991a (. Using the eight SPT DS¬ 


FGs with [Cii] and CO (1-0) detections, we determine the 
T[cii]/Tco(i-o) ra tio to be 5200 ± 1800. In section 


5.1 


derive the physical conditions that can be derived from this 
ratio from first principles and eventually compare it with 
the PDR models. 

Finally, in Figure|8j we test if the T[cii]/Tco(i-o) raticQ 
depends on Td. We find no correlation in either the low -2 or 
the SPT DFSG sample. However, the four high -2 sources 
with the lowest T[cn]/Tco(i-o) ratio are all AGN domi- 


9 Note that both parameters are independent of the lensing mag¬ 
nification. 
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Figure 7. The [Cn] luminosity versus the CO luminosity for the 
SPT DSFGs and the comparison low and high- 2 : samples. These 
star forming systems show a correlation between the [Cn] and 
CO (1-0) luminosities. Fitting a ratio to 11 SPT sources with [Cn] 
detections and CO (1-0) data yields a slope of ~ 5200. The width 
of the grey shaded area represents a la spread, ~ 5200 =b 1800. 
Fitting a ratio to the low- 2 : sample yields a slope of 1300 d= 440. 
The typical error bar for the low- 2 : sources is represented by the 
black cross. 



20 40 60 80 

Tdust [K] 


Figure 8. h[cil]/Tco(l 0 } vs Td f° r the SPT sources and the 
low and high -z comparison sample. The typical error bar is repre¬ 
sented by the black cross. The grey shaded area represents the I rr 
spread of the Tfcilj/hcofl-O) ratio in the SPT sample. Both the 
SPT sources and the low -2 sample do not show any dependence 
on Tu. The high -2 AGN-dominated sources are both warmer and 
have fainter [Cu] relative to CO. 


nated. This is consistent with the observations of |Stacey| 
|et al.| ( |2010[ > and |Sargsyan et al.| (2014| ) that AGN domi¬ 
nated sources have lower T[cii]/Tfir ratios. Although our 
Td determinations are rather crude, these AGN-dominated 
sources are also those with the warmest Td, as expected (see, 
e.g., Figure 3 in Greve et al. 20121. Interestingly, only 2 out 


of 11 SPT DSFGs fall near the AGN-dominated sources in 
Fig. i supporting our conclusion (§ |4.1.3| ) that the SPT sam¬ 
ple does not contain many strongly AGN-dominated sources. 


5 DISCUSSION 


5.1 Possible origins of the [Cii] to CO correlation 

5.1.1 Determining the optical depth of the [Cii] line 

The origin of the correlation between the [Cii] and CO (1-0) 
luminosities and the impact of their relative optical depths 


was first discussed by Crawford et al. (19851, and later by 

others in both low and high 

-2 objects (e.g. 

Wolfire et al. 

1989; Stacey et al. 1991a, 2010; |Swinbank et 

al.||2012; Neri 

et al. | 2014). 

Crawford et al. 

(1985) assumed [Cii] excita- 


ply the Rayleigh-Jeans approximation. This assumption was 
supported by independent estimates of the gas temperature 
of order ~ 300 K using other fine-structure lines and assum¬ 
ing optically thin [Cu] emission (Ellis fe Werner|1984| . They 
remarked that their observed T[cii]/^co(i-o) ratio ~ 4400 
is close to the [Cii] to CO(l-O) frequency ratio cubed, sug¬ 
gesting optically thick [Cu] emission. However, using sup¬ 
porting data in Orion, the Galactic centre and M82, they 
independetly derived [Cii] optical depths t[cii] = 0.03 — 1. 
This conclusion that [Cu] is mostly optically thin has been 


assumed several times since (e.g. 

Stacey et al. 1991a 

Hailey- 

Dunsheath et al. 2010 

Stacey et al. 2010; De Breuck et al. 

2011||Rawle et al. 2014 

)• 


The most accurate way to determine the optical depth 
of the [ 12 Cn] emission line is through observations of the iso¬ 
topic line ratios. The two [Cii] emission lines for 12 C and 13 C 
have been observed in local star forming regions, e.g. in M42 


(Stacey et al. 1991b Boreiko & Betz 19961, NGC2024 (in 
the Orion nebula; Graf et al.|2012 i, the Orion Bar, Mon R2, 
NGC3603, the Carina Nebula and NGC7023 (Ossenkopf 
et al. |2013 1. These observations find optical depths ranging 
from r ~ 1 to t ~ 3, with on average a moderate opti¬ 
cal dept t ~ 1.4. However, one should keep in mind that 
Orion and Carina are very bright [Cu] emitters (to allow for 
a detection of the faint 13 C), with very strong FUV fields 
(Go ~ 10 4 — 10 s ), implying large [Cii] columns and hence 
rather high opacities. The galaxy-scale average may result 
in lower optical depths near unity. Unfortunately, this iso¬ 
tope ratio technique cannot be applied to galaxy-integrated 
[Cii] observations as the [ 12 Cii] fine structure line and the 
brightest [ 13 Cii] hyperfine structure line are only separated 
by ~ llOkm/s, which is smaller than the typical widths of 
these lines. The much weaker 13 C line is then indistinguish¬ 
able from features in the profile of the 12 C line. 

The optical depth of [Cu] in very distant galaxies must 
therefore be estimated by other methods. The idea of opti¬ 
cally thick [Cii] emission at high redshifts was recently pro¬ 
posed by |Neri et al. (2014) for the high -2 sub-millimetre 
source HDF 850.1. Assuming the [Cii] line excitation tem¬ 
perature is the same as the dust and gas kinetic temper¬ 
atures, they argue for high line optical depth of the [Cii] 
line (i.e. rrcni > 1). The observed Lrcn] and predicted 


10 The [Cii] ground state energy level is 92 K. 
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ico(i-o) by |Walter et al.||2012| , yield a T[cii]/T C o(i-o) ra¬ 
tio of ~ 5200, in agreement with the ratio derived from the 
SPT sample of 5200 ± 1800 (§ |4.4| ). We now generalise this 
line ratio method by comparing the source functions. This 
does not require the CO and [Cii] excitation temperatures 
to be the same, nor that the Rayleigh-Jeans approximation 
applies (i.e. we also consider cases where T e x,[cu] is close to 
or below 91 K). 

From the source functions of both lines, the luminosity 
ratio depends on 
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where we have assumed that the [Cii] and CO emitting gas 
have the same filling factors (see j ]4.3| ). This assumption is 
consistent with the very similar [Cii] and CO velocity profiles 
(Figuretl I. In this case, Al i CI11 — _ _[ 4 cn]_ ^lso j n 

v |_| fll, CO(l-0) "COll-O) 

galaxies, [Cii] and CO have been found to trace each other 
both morphologically and kinematically (e.g. |Braine et al.| 
2012 and Mittal et al. 20111. Equation (2) has four free 
parameters: two excitation temperatures and two opacities. 
In order to obtain constraints, we therefore have to fix some 
of these. 


5.1.2 Same [Cii] and CO excitation temperatures 


We first consider the case suggested by Neri et al. (2014J) of 
equal excitation temperatures (i.e. T eXj [ C ii] = Tex,co(i-o)) 5 
Figure [9] plots T[cii]/^co(i-o) vs T ex for the three scenarios: 
i ) optically thin [Cii] and (nearly) optically thick CO(l- 
0) emission, ii ) (nearly) optically thick [Cii] and optically 
thin CO(l-O) emission, in) same optical depth for [Cii] and 
CO(l-O). 

Scenario i shown by the black dot-dashed curve under¬ 
predicts the observed ratio (grey shaded area) by an order of 
magnitude. In the opposite case, scenario ii (grey dot-dashed 
curve), the observed ratio is reached only for very low excita¬ 
tion temperatures. However, this optically thin CO scenario 
can be ruled out because both the 12 CO to 13 CO ratios and 
the low- J 12 CO line ratios in the SPT sample imply that CO 


is moderately optically thick with 7co(i-o) = 1 — 10 (Spilker 
etTal. 2014|), similar to what is seen in our own galaxyj fPen - 
zias et al.|1972|[Goldreich fc Kwan|1974| ). Scenario in, where 
the [Cii] and CO (1-0) optical depths are in the same regime 
is thus the only one that can fit the observed ratios, but 
only for excitation temperatures > 180 K. While we can¬ 
not distinguish mathematically between low and high opti¬ 
cal depth, the known rco(i-o)= 1—10 and the equality of 
both the excitation temperatures and opacities would imply 
that also [Cii] would need to be (nearly or fully) optically 
thick. 

The only way the above scenario in can fit the observed 
T[cn]/Tco(i-o) ratio is in a thermalised region where both 
[Cii] and CO(l-O) have excitation temperatures > 180 K 
(see Fig. |9j). However, the average CO excitation tempera¬ 
tures in the SPT sample are < 50K (Spilk er"et al.|20l4| , well 
below the values required to fit the observed ratio (Fig. |9]) . 
The only way to reconcile the model in equation (2) with 
the known physical parameters of the CO (tco(i-o) > I and 
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Figure 9. The T[cn]/Tco(l- 0 ) ratio predicted from the source 
functions versus equal [Cii] and CO(l—0) excitation temperatures 
(T ex , [C ii] = Tf.x,CO(i o) ) ■ f° r three different cases of the optical 
depth: i) optically thin [Cii] and optically thick CO(l—0) (the 
black dot-dashed curve), ii) optically thick [Cii] and optically thin 
CO(l-O) - (the grey dot-dashed curve), in) same optical depth 
of [Cii] and CO(l-O) (blue dashed curve). The grey shaded area 
represents the 1 rj spread of the T[cn]/Tco(i-0) ratio in the SPT 
sample. Case i underpredicts the ratio by an order of magnitude. 
Cases ii and in can both reproduce the observed ratio. However, 
we know from 12 CO to 13 CO ratios that CO is optically thick (e.g. 
|Spilker et al.|20 14), which rules out case ii. Only case Hi, implying 
optically thick CO and [Cii], is consistent with all observational 
data. 


T ex ,co )$ 50 K) is to allow for different excitation tempera¬ 
tures of [Cii] and CO. Such different excitation temperatures 
also imply different [Cii] and CO emitting regions within the 
GMCs. 


5.1.3 Different [Cii] and CO excitation temperatures 

To examine the model from equation (2) with differ¬ 
ent excitation temperatures, we have to fix the observed 
T[cn]/Tco(i-o) ratio and at least two parameters (two opti¬ 
cal depths or an optical depths and a temperature). We first 
consider the case when tjcii] = Tco(i-o), and plot T eX: co(i-o) 

(blue dashed line with 


10 


as a function of T eXj [cu] in Figure 
grey shaded area illustrating the observed T[cn]/Tco(i-o) 
range). We conclude that T ex> [cii] > T ex ,co(i-o) throughout. 
In the optically thick case, we can also allow for 1 ^ Ten] 
< t co(i—o) ■ The red curve in Figure [To| illustrates this for 
tco(i-o)=4. Any further increase of the difference between 
Rcii] and vcofi-o)! will also increase the difference in T ex . 
We warn that once both the [Cii] and CO become strongly 
optically thick (i.e. 1 <C tjcii] < 7co(i-o))> one can no longer 
determine any differences between the optical depths, and 
hence no longer determine T ex ren]- 

Alternatively, we can also fix the CO opacity and ex¬ 
citation temperature based on existing observations of the 
SPT DSFG sample, and determine which T eXi [cn] values are 
predicted for a given T[cii]- Assuming the molecular gas is 
traced by CO (1-0) and the dust is thermalised allows us 
to fix the T ex ,co(i-o) = T d ~ 35K ( Weifi et al.|[20 13j). This 
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100 




Figure 10. T ex co(i-0) 35 a function of T ex> [cni. I n cases, 
the excitation temperature of [Cii] is higher than for CO(l-O). 
The blue dashed curve with the light grey shaded area shows the 
observed T[cii]/Lco(i-0) = 5200 =t 1800 range of SPT DSFGs in 
the case of equal [Cii] and CO optical depths. The red continu¬ 
ous curve and hashed area illustrates that the difference between 
the excitation temperatures becomes even more significant for 
T[cn]=l and tco(i-o) = 4. The dot-dashed curve and the dark grey 
shaded area illustrate the case of T[cn]= 0.1 and tqo(i-o)— 1* 
The hatched horizontal area marks T e x,CO(i-0) = 35 =b 10K. 


Figure 11. The L[cii]/Lco(i-0) ra Lo as a function of the [Cii] 
excitation temperature for a fixed T ex co(i- 0 )—35 K. The ob¬ 
served Z/[cn] /^ 00 ( 1 - 0 ) ratio in the SPT DSFG sources (grey 
shaded area) is achieved within T eXj [cn] ~ 60 — 90 K for equal 
[Cii] and CO optical depth (blue curve). In the case of T[ch] 
= 1 and 700(1—0) — 4, the T[cii]/-Lco(i- 0 ) ra tio is achieved for 
T ex> [cii] — 35 — 110 K. The optically thin t"[cii] = 0.1 and nearly 
optically thick Tco(i-o) = 1 case is reached by [Cii] excitation 
temperatures in the range ~ 240 — 330K (black dashed curve). 


value is consistent with the T^ n determined from the stacked 


If 


ALMA spectrum of the SPT sample (Spilker et al. 2014). 
This T eX) co (l-o) =35 K case is illustrated in Figure 
T[ CII] = Tco(i-O): this would imply T ex> [ C n] ~ 60 - 90 K. 
Raising tco(i-o) = 4 while keeping tjcii] = 1 would im¬ 
ply r eXiI cii] ~ 80 - 110 K. The observed T[cn]/Tco(i-o) ra¬ 
tios can also be reproduced with optically thin [Cii] and 
(nearly) optically thick CO (t[cii] =0.1 and t CO (i-o) = 1) 
when the [Cii] excitation temperatures are ~ 240 — 330 K. 
A determination of T ex jen] is therefore needed to determine 
Tfciii ■ If the gas densities are higher than the critical density 
for [Cii] (2400 — 6100cm -3 for the above range of T ex> [cn]; 

20121, T eXj [ C ii]~ T(gas). The gas temper- 


Goldsmith et al. 


atures can be obtained from observed [Cii] / [OI] ratios (e.g. 
|Stacey et al.||1983| |Lord et al.||1996[ |Brauher et al.|[2008| ) , 
inferred through the peak [Cii] antenna temperatures (e.g. 


Graf et al. 2012 Ossenkopf et al. 20131, or by theoretical 
modelling (e.g. Kaufman ct al. 1999[). These studies obtain 


T( gas) = 100 — 500 K, which from Fig. 11 would imply [Cii 


optical depths ranging from 0.1 to 1. However, the densities 
we derive from a comparison with PDR models (see j ]5.2.1| ) 
are in the range 100 — 10 5 (TableJT with half of our sources 
below the [Cii] critical density, so the T e x,[cn]~ T(gas) may 
not be valid for a significant part of our sample. Hence, op¬ 
tical depths of 0.1 < r < 1 are consistent with the observed 
line ratios in the SPT DSFG sample. 

We also compare the model predictions in Figure [TT] 
with the low- and high- 2 : comparison samples. The lower av¬ 
erage T[cii]/ico(i-o) ratios in these samples (though with a 
much larger spread than for the SPT DSFGs) imply T e x,[cn] 
between ~ 30 and ~ 200 K. In particular for the low -2 sam¬ 
ple, optically thick [Cii] and CO emission would imply very 


low Tex, [cii] ~ 40 K, well below the ground state energy. The 
difference between the SPT DSFG and the low -2 comparison 
sample could therefore be ascribed to a lower optical depths 
in the nearby sources. In the context of PDR models, this 
can also be understood as a decrease of the Go /n ratio. Re¬ 
ducing this parameter implies smaller [Cii] emitting columns 
and a smaller effective optical depth. The lower Go/n ratio 
in the low -2 sample could be an effect of the lower far-UV 
fields found in galaxies forming stars at a more modest rate 
( |Stacey et al.|1991a||1993| [Kaufman et al.|1999| ). 

Finally, we note that Figure [8] also contains some low -2 
sources with T[cn]/Tco(i-o) ratios >15,000. Such values are 
difficult to explain with standard PDR models. Low metal- 


licity has been invoked to explain these sources (e.g. Mai 


oney fc Black|1988 Stacey et al.|1991a[|Madden et al.|1997 1. 
Considering this effect is beyond the scope of this paper. 

In summary, the observed T[cn]/Tco(i-o) ra ti° s i n the 
SPT DSFGs are best described by a non-uniform medium of 
[Cii] and CO(l-O) emitting gas with T e x,[cn]> T e x,co(i-o), 
Tco(i-o) 1 and tjcii] 1- 

5.2 Implications of different [Cii] and CO(l-O) 
excitation temperatures 


In section |5.1.2| we concluded that a homogeneous region 
with thermalised CO and [Cii] gas is incompatible with the 
observed T[cn]/Tco(i-o) ratios. Even in the optically thick 
case, the only way to reproduce the observed ratio is for 
higher [Cii] than CO(l-O) excitation temperatures (see Fig¬ 
ure [To]) . The cases of uniform and separated [Cii] and CO 
gas was recently studied by Mashian et al. ([2013). They ex¬ 
plore four different models to explain the observations in 
the high -2 submillimetre source HDF 850.1: 1) separate CO 
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- [Cii] virialised gas, 2) separate CO - [Cn] unvirialised gas, 
3) uniformly mixed CO - [Cii] virialised gas and 4) uni¬ 
formly mixed CO - [Cii] unvirialised gas. Based on cosmo¬ 
logical constraints due to the dark matter halo abundance 
in the standard A cold dark matter cosmology, they rule out 
three of the models and conclude that the preferred model 
is an unvirialised molecular cloud model with independent 
CO and [Cii] emitting gas with a average kinetic tempera¬ 
ture of 100 K and density of 10 3 cm -3 for the molecular gas. 


Both our conclusions and those of Mashian et al. (20131 are 


completely consistent with the structure described by PDR 
models. 

PDRs are clouds of molecular gas associated with star¬ 
forming regions, as they are often found near young massive 
O and B stars, acting as the source of the FUV photons that 
determine the temperature and chemical composition of the 
gas (Meijerink et al. 2007). Schematically, in a PDR the in¬ 
creasing extinction (Ay) with depth into the cloud creates a 
layer-structure, where the surface of the cloud with Av ~ 1 
is dominated by H + , C + and Oi. As the gas becomes more 
self-shielded against the dissociating FUV photons deeper 
in (Ay ~ 2 — 4), layers of Hi and H 2 form and a transi¬ 
tion region of C + , C and CO is present. At the centre of the 
cloud, the molecular gas is so opaque that the chemistry and 
heating are dominated by cosmic rays. Hence, the [Cii] fine 
structure line probes the surface of a PDR where Av < 1 
and T > 100 K, while CO traces the core of the cloud. To 
derive detailed physical parameters from the PDR models 
therefore requires spatially resolved observations of different 
species such as C, C + , CO, O, H 2 , polycyclic aromatic hy¬ 


drocarbons and dust continuum emission (e.g. Hollenbach 
fe Tielens|1999 Orr et ah|2014 1 . 


5.2.1 Implications of the photodissociation region 
structure 

The structure of the gas in PDRs allows for different [Cii] 
and CO(l-O) excitation temperatures. The physical param¬ 
eters for the ISM in low- and high -2 galaxies predicted by 
basic PDR models can be compared with observed data in 


a diagnostic diagram (Figure 12 I first presented by Wolfire 
et al. (19891 and updated by e.g. Stacey et al. (1991a 20101 
and |Hailey-Dunsheath et al.| ( |2010[ ). The diagram has the 
advantage of plotting ratios where the lensing magnifica¬ 
tion factors for high -2 sources, beam-filling factors for low -2 
sources and filling factors for the PDRs (which are unknown, 
but assumed the same for [Cii] and CO (1-0)) are divided 
out. It can be used to roughly estimate the strength of the 
FUV field and the gas density. One has to be cautious us¬ 
ing this diagram, as the FIR luminosity can contain strong 
contributions from other sources not associated with PDRs 
(e.g. AGN tori). 

The SPT sample is ideal in this respect as [Cii] has been 
measured for 17 out of 20, and low-J CO for 11 (see Table]!]) 
of these sources, so the CO luminosities for our sources (un¬ 
like the comparison samples) do not depend on uncertain 
scaling factors. The SPT sample is the most complete high- 
2 sample included in this diagram. The comparison sample 
is also integrated over entire galaxies enabling a fair com¬ 
parison to the SPT sample. 

The PDR model used in this diagram (Kaufman et al. 


zones of different depths. The intensities are modeled for a 
parameter space of the FUV field strength (Go) in units of 
the local Galactic interstellar ra diation held (the ‘Habing 
Field’, 1.6 x 10 -3 ergs cm _3 s _1 ; Kaufman et al. 19991 in 


the range 10~ ua ^ Go ^ 10°' 5 and the gas densities in the 
range 10cm^ 3 ^ n ^ 10' cm -3 . The emission from different 
species depends on the density, the field strength and the 


depth into the cloud. Hailey-Dunsheath et al. (2010) assume 


that 70% of the [Cii] emission originates from PDRs (see also 


5.2.21, meaning that the points would move slightly down 


in Figure [12| if this was corrected for. 

By comparing the SPT data points with the PDR 
model tracks in Figure |12[ we obtain a rough estimate 
of the radiation held strength and the gas density of 


100 < Go < 8000 and 10 2 cm“ 3 


< n < 10 5 cm 3 (see 


Table [2J. These values are consistent with the ones found 
in previous samples of DSFGs (e.g. Stacey et al. 20i0). 
They imply PDR surface temperatures of 300 — 500 K (Fig. 
2 of Kaufman et al. 19991. As these surface temperatures 
are representative for regions up to Av ~ 2, they cover 
most of the [Cii] emitting region. They are consistent with 
other derivations of the gas temperatures (see § 5.1.3). 


If the densitity exceeds the [Cii] critical density, these 
temperatures also represent T’exjcn]! which would imply 
optically thin [Cii] (Fig. 111. However, in half of our sources, 
the densities are below the critical densities, so a range of 
opacities up to unity remains possible. 

Sizes of PDRs 

Comparing the I/fir//i for the SPT source^ 33 ] which lie 
in the range (1.1 — 21.9)^t _1 x 10 12 Lq, with the local 
starburst galaxy M82 (Tfir~ (2.3 — 3.2) x 10 10 L 0 


et al. 1988 and Colbert et al. 19991 show the significant 


Rice 


difference in Lfir at high and low redshifts. Using this 
comparison and estimates of Go for each SPT DSFG, we 
estimate approximate sizes of the PDRs populating the 
galaxies, following |Stacey~et al.|pOlO). For this, we assume 
the molecular clouds are randomly mixed with young 
stellar clusters, acting as the radiation sources, within the 
galaxy (Wolh re et al.|[l990[ ). Assuming this structure, the 
relationship between the average Go, the total size of the 
PDRs ( D , diameter) and the total luminosity (Tfir) of 
the source is given by Go oc ALfir/D 3 for a short mean 
free path (A) and Go oc Lfir/D 2 for large mean free path 
of the FUV photons (see |Wolfire et aT. 19901. To estimate 
the approximate sizes, we read off Go and Lfir for the 
SPT DSFGs from Figure |12[ and scale these with the 
values obtained for M82 (Go ~ 1000 Lord et al. 1996 
2.8 x 1O 1O L 0 


Tfir^ 


average of the values determined by 


Rice et al. 1988 and Colbert et al. 19991, assuming the 


1999) models a plane-parallel slab divided into a number of 


same mean free path for the SPT sources and M82. The 
exact size of the PDR region in M82 is rather uncertain 
with reported sizes ranging from 300pc (Joy et al. : i 19871 to 
600pc (Carlstrom & Kronbcrg 19911. For consistency with 
|Stacey et al. (2010J), we will assume D ~ 300 pc. We warn 
that differential lensing (if significant, see § |3.2| ) could affect 
the positions of the sources in Figure [12] leading to different 
estimations of Go and therewith the sizes. The radii we 


11 See ( j .3.2 1 lor a discussion on the lensing magnification factor // 
in our sample. 

























































































14 


B. Gullberg et al. 



Lco(i-o/Lfir 


Wolfire 


Figure 12. T[ciii/Tfir vs T coCL-Cu/f'FlR. for the SPT sources and the low and high -2 comparison sample. The figure (e.g. 

|et al.|1989 , Stacey et al.|20lo| compares the values of the strength of the radiation field Go and the density n for low and high -2 sources. 
The diagram is independent of lensing magnification factors for high -2 source and beam filling factors for low -2 sources as both the [Cll] 
and CO(l-O) emission is normalised by the FIR luminosity. The typical error bar is represented by the black cross. To compare the 
observations with the model contours |Stacey et a l. (2010) assume that 70% of the [Cll] emission originate from PDRs. The grey shaded 
area represents the lcr spread of the Lr cn]/Tco(l-0) ratio in the SPT sample. 


estimate are listed in Table [2] along with our estimated Go 
and n. An estimate of the source radii from lens modelling 
is available for four sources in Table [2] The estimated radii 
of the PDRs, while fairly uncertain, are comparable to the 
source sizes determined by the lens models. 


A rough estimate of the size ranges occupied by the 
molecular gas is given in Table [2] as well. These sizes are 
estimated using the molecular gas mass estimated from 
Lqo(i-o)- Using low and high CO-to-fh conversion factors 
often used in the literature of aco = 0.8 Mg (K km/s pc 2 ) -1 
(see section |4.2| |Downes fe Solomon| (|1998 1) an d aco = 


Swinbank et al. 


2011), we esti- 


2M 0 (Kkm/spc““) _1 (e.g. 
mate sizes of the molecular gas regions making the simplistic 
assumption that the gas is uniformly distributed in a sphere 
with radius R and density n as listed in Table [2] The sizes 
estimated using this method are quite uncertain as we do 
not take into account non-uniform density profiles or non¬ 
unity volume filling factors of the gas. They are roughly con¬ 
sistent with both the sizes obtained from the PDR model, 
and the more accurate sizes derived from the lens models. 
We note that these kpc-scales are close to the typical sizes 
of host galaxies, and could be easily spatially resolved, es¬ 
pecially given the lensing magnification, allowing detailed 
future studies of their spatial distributions. 


5.2.2 Other contributions to the [Cn] emission 


The [Cii] emission integrated over an entire galaxy will con¬ 
tain contributions from regions with different physical con¬ 
ditions such as XDRs, CRDRs, shock dominated regions, 
diffuse warm gas, Hii regions, and PDRs. Above, we have 
assumed that observed [Cii] emission in the SPT DSFGs is 
dominated by emission from PDRs on molecular cloud sur¬ 
faces (e.g. Stutzki et al.|[l988 Stacey et al.||1993 l. We now 
consider the possible contributions from the alternative [Cii] 
emitting regions in increasing order of importance. 


In the vicinity of an AGN, supplying X-ray radiation, 
we expect XDRs. X-ray photons penetrate deeper into the 
volume of the interstellar clouds than the FUV photons in 
PDRs as the absorption cross sections are smaller for X- 
ray energies. As argued in § |4.1.3| and |4.4| the SPT sample 
does not show evidence for strong AGN activity, and XDR 
contributions are therefore expected to be negligible. 


In CRDRs, the gas heating and chemistry are controlled 
through interactions with high-energy particles. As the en¬ 
ergy density in cosmic rays is low compared to photons, 
CRDRs are thought to trace the dense, innermost regions of 
giant molecular clouds (GMCs), rather than the outer sur¬ 
faces where [Cii] emission is assumed to be more prevalent 
I Viti et al. 2013| ). We therefore assume the CRDR contribu¬ 
tion to be negligible. 

Turbulence and shocks have been suggested by | Applet] 
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Source 

Go 

n 

cm" 3 

R(PDR) 

kpc 

Ri /2 (le ns model) 
kpc 

^CO(l-O)) 

kpc 

h 

SPT0113-46 

10 3 

10 3 

0.6 - 

1.1 

— 

0.85 - 

1.16 

—t 

SPT0345-47 

10 4 

5 x 10 4 

0.4- 

0.7 

— 

0.39 - 

0.52 

—t 

SPT0346-52 

8 x 10 3 

10 5 

0.7- 

1.5 

0.59 + 0.03 

0.18 - 

0.24 

5.4+0.2 

SPT0243-49 

5 x 10 3 

10 3 

0.4- 

0.6 

— 

0.85 - 

1.15 

-t 

SPT0418-47 

10 3 

100 

0.7- 

1.5 

1.07 + 0.17 

2.26 - 

3.07 

21.0+3.5 

SPT0441-46 

8 x 10 3 

10 4 

0.3 - 

0.5 

— 

0.36 - 

0.49 

—t 

SPT2103-60 

10 3 

100 

0.7- 

1.4 

— 

1.69 - 

2.30 

—t 

SPT2146-55 

5 x 10 3 

10 4 

0.5 - 

0.8 

— 

0.42 - 

0.56 

—t 

SPT2147-50 

5 x 10 3 

10 4 

0.4- 

0.6 

— 

0.33 - 

0.45 

~t 

SPT0551-50 

10 2 

100 

1.0 - 

2.6 

— 

1.40 - 

1.89 

—t 

SPT0538-50 

10 2 

100 

1.7- 

5.9 

1.07 + 0.25 

2.46 - 

3.34 

20.9 + 4.2 


Table 2. The first column contains the names of the sources for which both [Cll] and low-J CO lines have been detected. The second 
and third columns list Go and n for the sources determined from the PDR models in Figure [I2| Note that especially n can be very 
uncertain as the models are very degenerate in this part of the diagram. Column four gives the size range determined using the short and 
long mean free path assumed by Stacey et al. (2010). The fifth column lists the radii for the sources which have lens models (Hezaveh 
et al. 2013), and the sixth column gives the range in sizes of the molecular gas estimated from t he molecular gas mas s range given by 
«CO = (0.8 — 2) Mq (Kkm/spc 2 ) -1 . The last column lists the lensing magnification factor from [Hezaveh e~ al. ( |2013[ ); SPT0529-54 is 
not included in this table as low-J CO lines have not been observed for this source. The sources marked with f, we have assumed a mean 
of ( fi ) = 14.1 (see j ]3.2^ . 


ton et al. (20131 to be an additional source of [Cn] emis¬ 
sion. They suggest that this mechanism should be present 
in highly turbulent conditions such as colliding galaxies and 
the early stages of galaxy-disc build up. However, for low 
turbulent velocities, it becomes difficult to distinguish be¬ 
tween this mechanism and low-density PDRs. The extreme 


L[cii]/£fir ratios Appleton et al. (2013) find for the inter- 
galactic filament in Stephan’s Quintet are ~ 30 x higher 
than those observed in the SPT sample. Smaller shock- 
ionised regions may still contribute significantly to the ob¬ 
served [Cii] emission. However, it would be difficult to ex¬ 
plain the relatively narrow spread in the observed [Cii] to 
CO ratio if a range of such shock-ionised regions would be a 
frequent occurence in our SPT sources. We therefore do not 
expect this mechanism to be important in our samples. 

The [Cii] emission could also originate from the diffuse 
warm low density medium in between the GMCs. Spatially 
resolved [Cii] and CO observations are required to differen¬ 
tiate between the PDR and diffuse components. Observa¬ 
tions of another lensed DSFGs, HLSJ091828.6+514223 at 


z = 5.243 (Rawle et al. 2014), do show that the [Cii] and 
CO(l-O) have a consistent structure and velocity profile, 
though the resolution may not go down to the scales of the 
GMCs (a few hundred pc). Observations of [Cii] and CO at 
spatial resolutions of ~50 mas are required to constrain this 
diffuse component in our DFSGs. 

Hu regions surround young O and B stars which emit 
Lyman continuum photons with energies exceeding the ion¬ 
isation energy of hydrogren (13.6eV). [Abel (20061 explore 
the contribution of [Cii] emission from Hu regions for a wide 
range in temperature, ionisation parameter ( U ) and electron 
density, and find that at least 10%, and sometimes up to 50- 
60% of the total [Cii] emission comes from within the Hii 
regions. Observations of other fine structure lines that only 
trace Hu regions are required to determine the exact con¬ 


tribution from Hii regions. Oberst et al. (2006 2011) used 


the observed [Cii]/[Nii] 205 /im ratio in the Carina nebula 
to constrain the contribution from Hii regions to 30%. This 


technique has since also been applied to high redshift objects 
( Ferkinhoff et al.|2011 Decarli et al.|20l4 i. 


5.3 The [Cii] /FIR luminosity deficit 

Several studies have reported a ‘deficit’ in the ratio of the 
[Cii] line strength to the FIR luminos ity ratio ( Ato m /Lfir 
for LIRGs with Lfir > 10 11 Lq (e.g. Malhotra et al.|[l997 


Luhman et al.|1998[|Maiolino et al.|2009[|Stacey et al.|2010 


Gracia-Carpio et al. 201 1[ Sargsyan et al. 2012| ). Various 
physical explanations for this trend have been proposed, in¬ 


cluding an increased ionisation parameter (Malhotra et al. 


2001 Abel et al.|2009| Gracia-Carpio et al.|2011 1, collisional 

de-excitation of [Cii] (Appleton et al. 20131 and non-PDR 
contributions to the FIR luminosity ( Luhman et al. ||2003[ ) , 
possibly from AGN (Sargsyan et al.|2012l. 


Stacey et al. (20101 argued that the lower [Cn]/FIR 


ratio can be explained by the fact that the star-formation 
in local ULIRGs is confined and vigorous (leading to high 
Go), while in the most distant objects, the star-formation 
is very large-scale, but of lower intensity (i.e. lower Go). In 
PDR models, the [Cn]/FIR ratio is inversely proportional 
with Go- Contrary to the FIR luminosity, which scales lin¬ 
early with Go, the [Cii] luminosity increases only slowly with 
Go- This is because in the observed density regime, the C + 
column density scales only with dust extinction, while the 
emissivity is only weakly dependent on Go since the gas 
temperature is above the excitation potential of 92 K. This 
leads to a “saturation effect” of the [Cii] emission at high 
luminosity in nearby ULIRGs, while the FIR remains unsat¬ 
urated. This is consistent with a recent study by |Dfaz-Santos| 
et al. (20141 exploring the difference in the [Cii]/FIR lumi¬ 
nosity deficit between the extended and compact (nuclei) 
regions in nearby LIRGs, revealing a larger deficit in the 
[Cii] /FIR luminosity ratio for the compact regions than for 
the extended regions. The ‘deficit’ is mostly confined to the 
innermost compact regions, while [Cii] to FIR luminosity 
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ratio for the extended regions is similar to that found in the 
extended disks of normal star-forming galaxies. 

A similar saturation effect would also occur when the 
[Cn] line becomes (nearly) optically thick. The line then 
reaches its maximum brightness, in the sense that any addi¬ 
tional incoming ionising photons will not increase the bright¬ 
ness of the line further. However the cooling of these addi¬ 
tional ionising photons may still continue through optically 
thin processes, notably the dust continuum. An alternative 
or additional way of decreasing the [Cn] to FIR luminosity 
ratio is therefore to have optically thick [Cii] and optically 
thin dust continuum emission. The higher density in the 
more compact regions may then increase the probability of 
the [Cii] to become optically thick. 

Finally we note that the [Cii] to FIR ‘deficit’ at high lu¬ 
minosities is also reported for other fine structure lines such 
as [Oi], [Om], [Nil] and [Nm], indicating that this ‘deficit’ is 
a general aspect of all FIR fine structure lines, regardless of 
their origin in the ionised or neutral phase of the interstellar 
medium (e.g. Gracia-Carpio et al.|20lf Farrah et al.|2013|. 


6 CONCLUSIONS 

We have presented the first uniformly selected [Cii] survey 
of lensed DSFGs covering the redshift range 2 = 2.1 — 5.7. 
We have detected [Cii] for 17 out of 20 sources, 11 of which 
are also observed and detected in low-J CO lines. This 
sample facilitates statistical studies of the ISM at high 
redshift. Our main results and conclusions are: 


1. We fit single or double Gaussian functions to the 
CO and [Cii] velocity profiles, and find consistent velocity 
profiles in 13 out of 14 CO detections with ALMA. This 
suggest that differential lensing is not significant in these 
cases, and is consistent with the idea that the [Cii] and 
CO(l-O) emitting gas are spatially associated. 


2. The line luminosity ratio of the [Cii] and CO (1-0) 
detections for the SPT sources is ~ 5200 ± 1800, which 
agrees with the first reported ratio by Crawford et al.| 
(19851. The values derived from the SPT sample are consis¬ 
tent with both low-* and high -2 comparison samples, but 
with significantly smaller dispersion. This is presumably 
due to the homogeneity of the SPT selection and followup 
observations and absence of any known AGN-dominated 
sources, which have lower [Cii] to CO(l-O) ratios. 


3. The SPT sample covers the same spread in the 
Li [cii] /Lfir ratio as the Lfir > 10 11 Lq sources in both 
our low -2 and high -2 comparison samples. AGN-dominated 
sources increase the scatter towards lower L[cii]/Afir ratios 
in the comparison sample. 


4. We investigate the origin of the [Cii] emission us¬ 
ing the observed L[cii]/^co(i-o) ratio, and conclude that 
the observed ratio is best described by a medium of [Cii] 
and CO(l-O) emitting gas with T eXi[CII ] > T e x,co(i-o), 
optically thick CO (tco(i-o) > 1), and low to moderate 
[Cii] optical depth (t[cii] < 1). The structure of PDRs 
allows for such different excitation temperatures of the [Cii] 
and CO(l-O) emitting gas. Interestingly the PDR models 


converge to this [Cii] to CO(l-O) ratio for densities below 
10 s cm -3 . 


5. We revisit the Z/[cii]/Lfir ‘deficit’ observed for 
sources with Lfir > 10 11 Lg, which has been explained 
as a “saturation effect” of the [Cii] emission in compact 
regions with higher Go factors. An alternative or additional 
explanation for this saturation effect is (nearly) optically 
thick [Cii] emission. In this case the [Cii] line becomes 
saturated and reaches the maximum [Cii] brightness, while 
cooling via the FIR continuum emission continues. The 
variation in the L[cii]/Lfir ratio is therefore dominated by 
the variation in Lfir rather than Lrcnl- 


6. We determine the FIR. luminosity for both the SPT 
sample and the comparison low and high -2 sample in a 
consistent way, adding 11 SPT DSFGs to the L[ C h]/Lfir 
vs L C o(i-o)/^fir plot in Figure 12 We compare the 


SPT sample with PDR. models and estimate the radiation 
field strength and average gas density to be in the range 
100 < Go < 8000 and 10 2 cm -3 < n < 10 s cm -3 . 


The reliability of [Cii] as a tracer of star formation 
rate (SFR) has been explored by e.g. De Looze et al. (2014). 
They show that [Cii] is a good tracer for the SFR except 
for low metallicity sources. Determining the metallicity 
using other fine structure lines such as [Nn], [On] and 
[Om], and the contribution to the [Cii] emission from Hii 
regions ( Croxall et al.|20 12l are therefore key in finding the 
most reliable tracer of SFR. in nearby and distant galaxies. 
Thus far, all high -2 fine structure line measurements have 
been unresolved, but thanks to lens shear of gravitationally 
lensed sources it will become possible to resolve structures 
down to 100 pc scales (Swinbank et al. 2010). Hence, for 
PDRs of the sizes estimated here, future observations with 
ALMA will be able to spatially resolve [Cii] emission and 
other fine structure lines and provide new insight into how 
[Cii] emission traces gas. 
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APPENDIX A: NEW REDSHIFTS 


All but two sources observed in [Cii] were published by Weifi 


et al. (20131, where the full source names are listed. The 


two new sources were selected from the list of Mocanu et al.l 
( |2013 ), and observed as part of our ongoing ALMA Cycle 1 
project to determine redshifts of additional SPT DSFGs. 
Table [An lists their full names and ALMA band 3 continuum 
positions. The redshift of SPT2319-55 is based on detections 


of CO(5-4) and CO(6-5), while the redshift of SPT2311-54 
is based on CO (5-4) confirmed by our APEX [Cii] detection 
(Figure [T]). 
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Source 


MCII] 
[10 9 L 0 ] 


L CO(1-0) 

[10® L 0 ] 


ipiR 
[10 12 L 0 ] 


T d 

[K] 


lensing 

magnification 


[Cn] and CO reference 


SMMJ2135-0102 (Eyelash) 2.33 5.8 ±1.3 0.8 ± 0.04 40.9 ±9.1 36.0 ± 2.3 32.5 ± 4.5 

SDP.130 2.63 <65 10.1 ±2.9 17.5 ± 3.6 38.6 ± 2.0 6±1 

SDP.81 3.04 275 ±39 22.0 ± 5.0 33.1 ± 6.3 40.0 ± 1.8 25 ± 7 

PSSJ2322±1944t 4.12 <8.0 6.1 ±2.6 28.6 ± 10.2 55.8 ±6.1 5.4 ±0.3 


Swinbank et al. 2010 Ivison et al. 2010 


Valtchanov et al. 2011 Frayer et al. 2011 


Valtchanov et al. 2011 


Benford et al. in prep., Carilli et al. 2002a 


BR1202-0725* 1 * 
HDF850.1 


4.69 

5.19 


15.8 ± 1.8 
7.5 ±0.8 


5.4 ±0.7 
2.2 ± 0.5 


34.1 ± 8.4 
1.7 ± 1.6 


55.1 ± 4.3 
29.9 ± 7.7 


Wagg et al.|2012 Carilli et al.|2002b 


1.5 - 1.7 


Frayer et al. 


2011 


SDP.141 

4.24 

60 ±9 

38.7 ± 4.6 

66.2 ± 17.7 

46.1 ± 3.2 

10 - 30 

Cox et al. 

2011 






BRI1335-0417t 

4.41 

15.7 ± 2.5 

4.4 ±0.8 

17.7 ± 4.9 

51.3 ±4.4 

— 

Wagg et a 

1.1201 

0 

Cai 

illi et al. 

2002b 


ALESS65.1 

4.45 

3.2 ±0.4 

< 1.11 

3.12 ± 1.40 

44.3 ± 5.1 

— 

Swinbank et al. 


012 

Huynh e 

al.fet 

)14 


Neri et al. 2014 Walter et al. 2012 


HLSJ091828.6±514223 

5.24 

82.4 ± 2.3 

43.0 ± 3.3 

97.5 ±24.4 

48.7 ± 3.0 

8.9 ± 1.9 

Rawle et al. 

2014 







SDSSJ1044-0125t 

5.78 

1.6 ±0.4 

1.6 ±0.3 

5.6 ±5.1 

59.8 ± 25.8 

— 

Wang et al. 

2013 







SDSSJ2310±1855t 

6.00 

8.7 ± 1.4 

12.1 ± 1.0 

20.0 ± 14.5 

61.8 ± 20.7 

— 

Wang et al. 

2013 







HerMESFLS3 

6.34 

14.9 ±3.1 

4.6 ± 1.5 

29.9 ±9.6 

52.8 ±6.2 

2.2 ±0.3 

Riechers et 

il. 20 

13 1 






SDSSJ1148±5251 t 

6.42 

4.6 ±0.5 

6.4 ±0.7 

9.8 ±6.7 

60.9 ± 15.7 

— 

Maiolino et 

al. 20 

05 






IRASF10026±4949 

1.12 

25.8 ±3.2 

< 0.8 

— 

— 


Stacey et al. 

201C 







3C368 

1.13 

9.3 ± 1.5 

< 1.9 

— 

— 


Stacey et al. 

201C 

Evans et al.|l996 



SMMJ123634.51±621241.0 

1.22 

14.4 ± 2.0 

3.76 ±0.12 

— 

— 


Stacey et al. 

201C 

Frayer et al."|2008 



MIPSJ142824.0±352619 

1.32 

10.5 ±3.1 

1.10 ±0.03 

— 

— 


Hailey-Duns 

ieatb 

et 

al. 2010 

Iono e 

al.L 

006 

SDSSJ100038.01±020822.4t 

1.83 

10.6 ± 2.4 

3.10 ±0.17 

— 

— 


Stacey et al. 

201C 

Aravena e 

al. 

2008 


SWIREJ104704.97±592332.3 

1.95 

20.0 ± 2.6 

8.9 

— 

— 


Stacey et al. 

201C 







SWIRE J104738.32±591010.0 

1.96 

12.1 ± 3.2 

3.5 

— 

— 


Stacey et al. 

201C 







BRI0952-0115 t 

4.44 

4.5 ±2.6 

0.43 ±0.11 

— 

— 


Maiolino et 

al.|20 

09 

Guilloteau et al.|l99 

9 

LESSJ033229.4-275619+ 

4.76 

9.9 ± 1.5 

1.05 ± 0.25 

— 

— 


De Breuck e 

al.j 

2011 

Coppin et al. 

2010 


SDSSJ0129-0035t 

5.78 

1.8 ±0.3 

2.8 ±0.5 

— 

— 


Wang et al. 

2013 








SDSSJ2054-0005t 

6.04 

3.3 ±0.5 

2.7 ± 0.6 

— 

— 


Wang et al. 

2013 







ULASJ1319±0950t 

6.13 

4.4 ± 0.9 

3.5 ± 0.7 

— 

— 


Wang et al. 

2013 







CFHQSJ0210-0456 t 

6.43 

0.30 ±0.04 

< 0.88 

— 

— 


Willott et a] 

J 201 

3 






ALESS61.1 

4.42 

1.48 ±0.23 

— 

3.0 ± 1.4 

43.7 ± 5.1 

— 

Swinbank et 

al. 2 

012 



Table Bl. The high- 2 : sources in the comparison sample. All but one sources (ALESS61.1) have published [Cn] and CO detections. 
The observed CO luminosities have been scaled to CO(l-O) luminosities using the ratios from Stacey et al. 2010. The first 14 sources 
have a sufficient amount of photometric data published for the determination of Lp ir and X^, while the remaining are unconstrained. 
ALESS61.1 has good photometry, but no published low-J CO observations. 

* A sum of the [CII] and CO emission from the north and south source. 

' AGN dominated source. 


Short name 

Source 

R.A. 

Dec. 



J2000 

SPT2311-54 

SPT-S J231124-5450.5 

23:11:23.94 

-54:50:30.0 

SPT2319-55 

SPT-S J231922-5557.9 

23:19:21.67 

-55:57:57.8 


Table Al. Source names are based on positions measured with 
the SPT (Mocanu et al. 2013). Source positions are based on the 
ALMA 3 mm continuum data. 


APPENDIX B: HIGH-Z COMPARISON 
GALAXIES 


This paper has been typeset from a T^)K/ DT^gX file prepared 
by the author. 


























































































































